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ABSTRACT 



The strain response of metal plates during welding is 
discussed and current state-of-the-art efforts to analyze 
the phenomena are reviewed. Mechanical and physical temp- 
erature dependent properties; transient strain and temper- 
ature distribution data obtained from experiments on 
6061 aluminum alloy in the T651 condition; and a data 
reduction computer program are presented. The experiments 
were designed to look at the macroscopic effects of welding 
upon heat treatable age hardened alloys while at the same 
time approximating ship structural weldments. 

The transient strain response was found to be pre- 
dominantly longitudinal but transverse strains were signifi- 
cant in the region of the welding arc. It was also found 
that residual strains correlated reasonably between similar 
experiments with a variation of the second stress deviator 
tensor invariant J 2 representing the "state of stress" 
rather than by individual strain components. 

Several recommendations are made concerning continued 
experimental investigation aimed at further development 
of the National Aeronautics and Space Adminstration programs. 
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I: INTRODUCTION 



A. Background : 

The phenomenon of thermal stresses and strains resulting 
from welding processes lead to difficulties in the fabrication 
of structures. Non-uniform temperature distribution leads 
to residual stresses which in turn provide problems from 
stress corrosion, buckling, brittle fracture, fatigue strength, 
and distortion. Depending upon the type of material being 
welded, any or all of the above problems may be important 
considerations in the structure's design and its service per- 
formance (1) . All contribute to the reliability of the 
basic structure (2, pp. 2). Masubuchi (2 - 5), Klein and 
Masubuchi (6), and Klein (7) provide interpretive reports 
which review the particulars of each of the above phenomena 
for aluminum and steel as well as comprehensive reviews of 
the literature. Almost all of the previous studies refer 
to residual stresses and distortion and only recently has 
much attention been given to the actual thermal stresses 
during welding. This is due to the complexity of the problem 
characterized by large temperature changes in small areas 
near the welding arc with its resulting non-elastic deforma- 
tion, temperature dependency of the material properties and/or 
phase transformations, and complex boundary conditions 
resulting from conditions of geometry and multi-pass welding 
(3, pp. 2, 3) . 



12 

B. Thermal Stresses During Welding , Residual Stresses and 
Distortion : 

In order to study the effects of residual stresses 
and distortion performance of welded structures, it is first 
necessary to understand the relative magnitude and distribution 
of residual stresses and distortion in weldments. The 
mechanism by which thermal stresses and strains are developed 
in welded plates is best described by Masubuchi (8,3, pp. 4-6) 
and is repeated here: 

Figure 1 shows schematically changes of temperature 
and stresses during welding. A bead-on-plate weld is being 
made along the x-axis. The welding arc, which is moving at 
a speed, v, is presently located at the origin, 0, as shown 
in Figure la. 

Figure lb shows temperature distribution along several 
cross section. Along Section A-A, which is ahead of the 
welding arc, the temperature change due to welding, AT, is 
almost zero. Along Section B-B, which crosses the welding 
arc, the temperature distribution is very steep. Along 
Section C-C, which is some distance behind the welding arc, 
the distribution of temperature change is as shown in 
Figure lb-3. Along Section D-D, which is very far from the 
welding arc, the temperature change due to welding again 
diminishes . 

Figure lc shows the distribution of stresses along 
these sections in the x-direction, a . Stress in the 

X 
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y-direction, a , and shearing stress, t^, also exists 



in a two-dimensional stress field.* 

Along Section A-A, thermal stresses due to welding are 
almost zero. The stress distribution along Section B-B is 
shown in Figure lc-2. Stresses in regions underneath the 
welding arc are close to zero, because molten metal does 
not support loads. Stresses in regions somewhat away from 
the arc are compressive, because the expansion of these 
areas is restrained by surrounding metal that is at lower 
temperatures. Since the temperatures of these areas are 
quite high and the yield strength of the material is low, 
stresses in these areas are as high as the yield strength 
of material at corresponding temperatures. The magnitude 
of compressive stress passes through a maximum with increas- 
ing distance from the weld or with decreasing temperature. 
However, stresses in areas away from the weld are tensile 
and balance with compressive stresses in areas near the 
weld. In other words, 



across Section B-B.** Thus, the stress distribution along 



fa 



x 



dY = 0 



(1) 



* 




In general three-dimensional stress field, six stress 

components, rr„, o > a > t , t , t exist. 

u x u y u z l xy 1 zy l zx 



Equation (1) neglects the effect of a 
equilibrium condition. 



y 



T 



xy 



and 



on the 



14 



Section B-B is as shown in Figure lc-2. 

Stresses are distributed along Section C-C as shown in 
Figure lc-3. Since the weld metal and base metal regions 
near the weld have cooled, they try to shrink causing 
tensile stresses in regions close to the weld. As the 
distance from the weld increases, the stresses first change 
to compressive and then become tensile. 

Figure lc-4 shows the stress distribution along Section 
D-D. High tensile stresses are produced in regions near 
the weld, while compressive stresses are produced in regions 
away from the weld. The distribution of residual stresses 
that remain after welding is completed as is shown in the 
figure. 

The cross-hatched area, M-M, in Figure la shows the 
region where plastic deformation occurs during the welding 
thermal cycle. The ellipse near the origin 0 indicates the 
region where the metal is melted. The region outside the 
cross-hatched area remains elastic during the entire weld- 
ing thermal cycle. 

As shown in Figure 1, thermal stresses during welding 
are produced by a complex mechanism which involves plastic 
deformation at a wide range of temperatures from room temp- 
erature up to the melting temperature. Because of the 
difficulty in analyzing plastic deformation, especially at 
elevated temperatures, mathematical analyses were limited 
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b. Tamparotur* Chong* c. Strata <r % 



Figura l Schamatic Rapraaontatton of Chongaa of Ttonporotur* and Stroaaos 
Ouring WatdWtg 0 . 8 ). 
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for very simple cases such as spot welding. 

However, on the basis of recent developments in com- 
puter technology, it appears that a technological break- 
through as far as the analysis of thermal stresses during 
welding are concerned is not too far away. 

C. Previous Investigations : 

The closest approach to a breakthrough to date has 
been the one-dimensional strip theory approach first pre- 
sented by Tall (9) and extended by Masubuchi and associates 
at Battelle Memorial Institute (10) and Massachusetts 
Institute of Technology (7, 11, 12) under the support of 
NASA's George C. Marshall Space Flight Center. These 
studies have been reviewed in detail in references (3- 6) 
and will not be repeated here. 

A great deal of effort in the last several years has 
been devoted to developing finite element programs which 
will extend one-dimensional analysis to two and three 
dimensions (13 - 20) . Nearly all have met some qualitative 
success but in general, quantitative results have been 
lacking. The finite element method utilizes an approximate 
model to derive a set of equations which is then solved 
exactly. The flexibility of the finite element method will 
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allow the extension of the flat plate problem into more 
complex and real life geometries. It provides a complete 
strain picture, including longitudinal, transverse, and 
shear strains. It has the capability to model the restraint 
conditions in the actual structure and weld. Its greatest 
drawback is high expense from considerable computer time 
for large degree of freedom problems and the post analysis 
of the vast quantities of data the technique generates. 
Interactive graphics will greatly aid this latter problem 
in the future. 

Masubuchi and Iwaki (13) and Masubuchi and Andrews (14) 
coupled the elastic-plastic analysis developed by Wang (15) 
with the thermal loading calculations for moving heat 
sources previously utilized by the one-dimensional programs. 
Reference (13) was qualitative and reference (14) provided 
some quantitative results with course meshes but appeared to 
have transition element instability when a smaller mesh was 
utilized. In this situation, an element which reaches 
yield stress is considered a transition element. However, 
in the next iteration, the equivalent stress from that 
element drops below yield and the element becomes elastic 
again. The whole situation jumps between the two sets of 
values. In order to overcome this difficulty, a scheme 
was devised to quarter the time step in which this insta- 
bility occurred. This worked for a large number of cases 
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but not all cases. Sometimes increasing the time step 
also eliminated the instability. Reference (10) may offer 
a way to correct this problem. Yamada (17) seems to have 
met with more success dealing with thermal stresses due to 
rapid heating by controlling the number of elements which 
yield. 

Hibbit and Marcal (18) offer the most analytic of the 
approaches to date. Their model treats the weld process 
as a thermo-mechanical problem. A finite element formula- 
tion derived from the uncoupled thermal and mechanical 
energy balances forms the basis of the model. It attempts 
to deal with material non-linearity due to temperature 
dependence of thermal properties and in the fusion problem 
where material phase change is accompanied by a latent heat 
effect. It includes radiation as a cooldown mechanism and 
finite strain effects during isothermal loading. However, 
there was little agreement with experimentally measured 
residual stresses and they concluded that their finite 
element model did not include significant material behavior. 

Shinn (19) approaches the two-dimensional distortion 
of a panel structure due to welding with the assumption 
of elastic deformation during the welding process. Com- 
putations were carried out using the one-dimensional 
experimental values of unconstrained angular changes along 
the welded edge and its equivalent constrained welding 
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moment as an input to the computer program. The results 
were not completely successful because of the elastic 
assumption and the questionable accuracy of the experimental 
data input into the program. It is of the opinion of this 
author that this technique, on the basis of current tech- 
nology, will provide the most fruitful results in the 
immediate future in terms of financial and manpower invest- 
ments. The basic welding problem is highly non-linear 
and may be too complex to solve in a completely analytical 
form. A purely empirical approach does not appear very 
fruitful, either. It is not a simple task to determine 
thermal stresses in small regions heated to high temperatures. 
Without proper analysis, it would be difficult to adequately 
interpret experimental data. 

The empirical-analytic technique, however, is only as 
good as experimental data provided. At present, material 
and physical properties at the range of temperature from 
room to melting, even for the most common of structural 
materials, are sorely lacking in the literature. More 
effort must be expended in providing this critical information. 
The approximations utilized in the interim will be the 
controlling feature of this computational technique. 
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D. Aim and Purpose of the Present Study : 

Muraki and Masubuchi are presently developing a new 
elasto-plastic finite element computer program at Massachusetts 
Institute of Technology along the lines described by Yamada 
(17) with temperature distribution adopted from reference (20) . 
It is the purpose of this report to provide experimental 
data on a heat treatable tempered aluminum alloy to verify 
this program. Electric resistance thermocouples, three 
element strain gage rosetts, single element strain gages, 
and an extensioneter will be utilized to obtain transient 
and residual strain and displacement data in the center and 
near the edge of 30x18 inch panels for an automatic GMA 
welding process. Experiments will be performed for bead-on- 
plate and butt welding processes under similar heat inputs 
to verify: 

1. that experiments may be duplicated under reasonably 
similar conditions to give repeatable results; 

2. to compare butt and bead-on-plate welds for similar- 
ities and differences; 

3. to compare strains in both top and bottom of the 
plates to study the importance of bending strains as com- 
pared with transverse and longitudinal strains; 

4. to observe transient principal strain magnitudes 
and directions; 

5. to measure residual stress/strain distributions 
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in the longitudinal and transverse directions; 

6. to measure the deflection of the two half plates 
in a butt welding process under tack welded conditions; 

7. to present longitudinal and transverse strains for 
computational technique verification, and; 

8. to collect and present in usable form available 
physical and mechanical temperature dependent data for 
use in the computer analysis. 

Standard strain gages will be utilized which may not 
exceed 400°F. Since the material properties in the heat 
affected zone (HAZ) are not fully known or understood and 
the high temperature strain gages have proven unreliable 
(6, 7) in part because of temperature compensation data and 
in part by the lack of material properties data within the 
HAZ, strain measurements will be beyond the HAZ. 
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II: MATERIAL BEHAVIOR 



A. General : 

A new finite element program to study residual stresses 
and metal movement is presently being developed by M.I.T.'s 
Department of Ocean Engineering concurrent with this thesis. 
Since an attempt is being made here to provide material 
properties for and consistent with this new program, a 
brief and simple formulation of the tangent stiffness 
technique, incremental stress strain relationship, and 
thermal loading terms is provided for the sole purpose of 
demonstrating the importance of temperature dependency. Then, 
material properties consistent with this formulation are 
presented for 6061 aluminum alloy in the T6 and T651 condi- 
tion, the material utilized for the experiments presented 
in Chapter III. 

B. Mathematical Formulation: 

The most common finite element approach to plasticity 
problems in welding has been the Tangent Stiffness Approach. 
The problem is first solved elastically for the general 
loads. Then the elastic limit loads are determined by 
scaling down the given loads until the element with the 
maximum equivalent stress is just reached. The differences 
between the given loads and the elastic limit loads are then 
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£ 

Figure 2: Equivalent Stress and Strain (14, pp 20) 




Figure 3: Iteration Process to Bring Back the Stress Point 

on to the <T- £ curve with H*. The element in the tran- 
sition region starts out at point B as elastic. The elastic 
increment BC of the element for the prescribed load increment 
dP is first calculated. The stress point corresponding to 
C should be D on the stress-strain curve. Using the modulus 
corresponding to BD, the increment BF for dP is next calcu- 
lated. This process may be repeated until the solution 
converges (17, pp 300,303). 
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divided into a predetermined number of loading increments. 
The problem is then solved including plasticity for each 
element. When an element becomes plastic, the tangent 
stiffness approach is used to determine an effective elasto- 
plastic stiffness for that element. In this approach, 
the Prandtl-Reuss relation and the plastic modulus are 
combined with the elastic stress-strain relations to produce 
an incremental stress-strain relation which includes the 
effects of plasticity. These incremental relations depend 
upon the instantaneous stress levels in each element (14, 
pp. 4 - 5) . 

In general, 

do = E, de (2) 



where E is the tangent modulus coefficient. E, has an 

✓sir - 



elastic term and a plastic term. The plastic term is zero 



if no yielding takes place. Then, E fc may be used to estab- 



lish the tangent stiffness array K for each element, 

'>5 




[ bJ; E. B d 

J , oB ^ ^n 

Volume 



(Volume) 



(3) 



where B defines the strain displacement relationships. 

Thus, where loading is applied in increments, the structure 
stiffness array may be found at each stage of loading by 
evaluating each element's tangent stiffness (21, pp. 12 - 13). 
The incremental pseudo-thermal loading term likewise 



becomes 



dP 



n 



■ ( 



B*" E. de d (Volume) 
Volume 



(4) 



where d^ is the thermal strain caused by the welding 
process. Summed over the entire structure, one obtains 



K dU 
dU = 
de = 
da = 



= ^dP 
jT 1 dP 
B t dU 



E fc B dU 



(5) 

( 6 ) 

(7) 

(8) 



To determine the tangent modulus coefficient related 
to equation (2) , one must consider the following relations 

Prandtl-Reuss Flow Rule: 



, p 8a ,-p 

= 30 de 

Yielding Surface Behavior: 

d3 ■ If 

Material Behavior: 



(9) 



( 10 ) 



o = o, 



+ Ep IT = o. 



da = E? de. 



t -p 
= H' de. 



+ H' e_ 



(ID 

( 12 ) 



Utilizing equations (9) - (12) and recalling that 



The p and e superscripts refer to plastic and elastic 
respectively. 



da = Dde e = D (de - de P ) ( 13 ) 

^ ''"S-*- **** ^s/ /^s_/ 

where de is the total incremental strain and D is given in 

'•w 

Equation (24) for the plane stress situation, one arrives 
at the following relationship: 
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da = 



D - D 




3a 

3a 


D 


m 


1 


JL £ 3a 

/%✓ 


. h1 


+ 


M t d 

8a ~ 


• 

1 lolo} 

<ro |<TD 


► 

3 



► de 



(14) 



This reduces to equation (2) where the quantity in ^ 
is the tangent modulus coefficient. This is consistent 
with reference (17, pp. 310) . 

If we assume plane stress conditions and consider the 
equivalent stress for Von Mises' isotropic material. 



-2 2 2,^2 

a = a - aa + a +3 t 

x x y y xy 



(15) 



Differentiating, 
da 



= (2a - a ) da + (2a 

2 - L x y' x v y 



This may be rewritten in matrix form as 



IT 

da = — a 



1 0 

’I 1 0 

0 0 3 



da 



or 



_ IT 

da = — a 1 F da 



a ) da + 6 t dt 
x ; y xy xy] 

(16) 



(17) 



( 18 ) 
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where 



a 



«_ 




— 


a 




da 


X 




X 


a_. 


(19) da = 


da„ 


Y 






T 




dx 


*Y 

* 




xy 



(20) 



Comparing equations (10) and (18) , we see that 

r-K' 

Transposing, 



( 21 ) 



3a _ 1 _T _ 1 _ 

__ = —Fa = — Fa 

da ~ ^ 



a a 



( 22 ) 



Substituting the previous relationships into equation (14) 
we arrive at 



da 



D - 



D F a a T F D 

A/ /y/ A/ A/ ^ 

o 2 h' ” 



+ a F D F a 

A/ /V Ay 



de 



(23) 



for a two dimensional, plane stress, isotropic case, with 

the Von Mises' yield criterion where 

1 
0 



D = 



V 

1 



0 

0 



F = 



1 - V‘ 



* -i 

-i 1 
0 0 



0 0 

0 

0 
3 



1-v 

2 



(24) 



-2 2 . „ 2 , - 2 
a = a x " a x a y + a y + 3T xy 



(25) 



( 15 ) 
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Equation (23) is but a simple example of the more 
complex general case. The important point to note here is 
that E, o, and H' are functions of temperature and this 
makes quite complex even the most simple of welding problem 
examples. 

In the previous study at the Massachusetts Institute 
of Technology, Andrews (14, pp. 25) assumed the simplest 
model possible for the nodal force matrix, i.e., that of 
a constant temperature of each element. In this case, the 
nodal force matrix for each increment becomes 



where 



dP 



n 



■ I 



(B E de ) d (Volume) 



Volume 



de. 



= AT 



a 

a 

0 



(4) 



AT = average temperature change of element 
a = coefficient of linear thermal expansion 



If the element in question is plastic, the tangent stiffness 
matrix is used. Masubuchi and Iwaki (13, pp. 9) assumed 
a much more complex loading function which accounted for 
the equivalent nodal forces due to the dependency of the 
coefficient of thermal expansion and Young's modulus on 
temperature. The plastic terms additionally took the 
temperature dependency of the yield surface into account. 
Yielding is complex. It can occur during the heating stage 
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and also in the subsequent cooling stage since some 
elements can be stressed to their yield point in the 
reverse direction. Yamada (17, pp. 298) assumed the same 
simplified type of expression as equation (4) in his 
study of thermal stresses due to rapid heating. As 
in the case of the incremental stress-strain relation, the 
important point to note is that in varying degrees of 
difficulty, the loading term is dependent upon E, a, and 
a which are in turn dependent upon temperature. 

The temperature distribution, uncoupled from the 
mechanical problem just described, may be approximated by 
classical, finite difference, or finite element techniques 
(9-11, 13, 14, 18, 20, 24-30). In terms of the conductiv- 
ity tensor, k, the general representation of the aniso- 
tropic heterogeneous continuum is 



where p is the density, c is the specific heat, k is the 
thermal conductivity, and U 1 " is the heat generation per 
unit volume (23, pp. 29, 44-45). Fourier's law applied 
to the boundary conditions provides 




V • (k • VT) + U'" 



(25) 



q 



— K • VT 



(26) 
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where q is the heat flux normal to the boundary surface. 
If the continuum is isotropic, these vector equations 
reduce to the familiar 

dT 

pc = V • (kVT) + U' M (27) 

q = -kVT (28) 

Equation (27) may be reduced to 



pc 



dT 

dt 



Vk 



VT 



kV 2 T 



U 1 



(29) 



If k is a function of space only, equation (29) is 
linear. On the other hand, where < depends on temperature 



alone , 



equation (29) may be modified to: 



dT 
PC dt 



dU 

dT 



(VT) 



kV 2 T 



U' 



(30) 



which is nonlinear. For homogeneous isotropic continua k 
is constant and equation (29) reduces to 



dT 

dt 



2 

aV T 



+ 



U m 

pc 



(31) 



where a is the thermal diffusivity. Since p, c, and k 
are functions of temperature, the welding problem becomes 
highly nonlinear. Coupled with phase changes, this 
problem becomes even more complex. 



C. Material Properties: 

As discussed in Section B, tabulation of the material 
properties as a function of temperature is of prime impor- 
tance. This information is difficult to obtain and often 
simply is not available. In this experimental work, alloy 
6061 in the T651 Condition is utilized. It is one of the 
most versatile of the wrought, heat treatable aluminum 
alloys. An understanding of the behavior of this common, 
relatively inexpensive material during welding and heat 
treating may also contribute to the understanding of the 
behavior of quenched and tempered steels such as HY-80 and 
HY-130 currently in use today. 6061 is readily welded by 
all methods and has excellent weldability characteristics. 



The plate used in this experiment has a Federal Specifi- 
cation of QQA-2 50/11. 

1. Composition (31, Code 3206, pp. 1) 
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Percent 
Min. Max. 


Chromium 


0.15 


0.35 


Copper 


0.15 


O 

• 

o 


Iron 


- 


0.7 


Magnesium 


0.8 


1.2 


Manganese 


- 


0.15 


Silicon 


0.40 


0.8 


Titanium 


- 


0.15 


Zinc 


- 


0.25 


Other impurities 
each 


__ 


0.05 


total 


- 


0.15 


Aluminum 


Balance 



TABLE 1: Composition of 6061 Aluminum Alloy 

2. Heat Treatment (31, Code 3206, pp. 1) 

Annealed (0 Condition) : Heat to 775°F for 2 to 6 

hours, cool at 50°F per hour maximum to 500 °F. The rate 
of subsequent cooling is unimportant. 

T4 Condition : Solution heat treat to 970°F, water 

quench, and naturally age at room temperature to a sub- 
stantially stable condition. 

T6 Condition : Solution heat treat to 970 °F, water 

quench, artificially age by precipitation heat treatment 
at 320°F for 16 to 20 hours or 350°F 6 to 10 hours and 
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air cool. 

T651/T451 Condition : Same as the T6/T4 Condition followed 

by a 1 1/2 to 3 percent permanent set stretch stress relief 
prior to any precipitation heat treatment. 

The soft, as quenched 0 Condition can be preserved 
by refrigeration in order to minimize springback and 
increase the general ease of forming operations. Maximum 
holding times which will preserve the formability of this 
condition are shown in Table 2 (31, Code 3206, pp. 3). 



Temperature °F 


R.T. 


32 


20 


Time 


2 hours 


2 days 


7 days + 



TABLE 2: Holding Times for 6061 in 0 Condition 



3. Specified Mechanical Properties (31, Code 3206, pp. 2) 



TABLE 3: Specified Mechanical Properties. 



Condition 


0 


T4, T451 


T6, T651 


Thickness - in. 


0.250 


0.250 


0.250 


Ultimate Stress in 
Tension min. ksi 


22 


30 


42 


0.2% Offset Yield 
Stress in Tension 
min - ksi 


12 


16 


35 


Ultimate Elongation 
2 inch min - % 


18 


16 - 18 


6-10 


Ultimate Shear Stress 
Typical - ksi 


12 


24 


30 
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4 . Mechanical Properties : 

In welding aluminum, it must be remembered that 
the metal in the area of the weld will be in an annealed 
condition after welding. There will be a corresponding 
loss of strength unless heat treatment is done after 
welding. Welding does not always reduce the strength 
of the heat treated alloys to that of the fully annealed 
condition because there is some air quenching of the metal 
as it cools (32, pp. 6.33). Table 4 (32/ pp. 6.34) 
demonstrates this effect for the type of filler metal 
utilized in this experiment. 



Filler 

Wire 


Specified minimum 
tensile strength of 
base plate in ksi. 


Average tensile 
strength across 
weld in ksi. 


Average bend 
free elonga- 
tion % 




Not heat treated after welding 


4043 


42.0 


27.2 


16.0 




Heat treated and aged after welding 


4043 


42.0 


43.5 


11.0 



TABLE 4: Strength and Ductility of Welded Butt Joints in 

Aluminum (TIG and GMA with argon) . 



Brungraber and Nelson (33) report that 6061 alloy may 
be heated up to 550° F by butt welding without appreciably 
affecting the mechanical properties of the base plate. 

They further report that for thicknesses less than one 
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inch, the distance from the center of the weld to the 
edge of the heat affected zone (HAZ) may be predicted 
as : 

b h = 0 . 6t + 0.00059 (32) 

where E has units of volts, I amperes, V inches per 

minute, and t (thickness) inches. As an example, for 

E = 20., I = 240., t = .25 and V = 32.16, El/Vt = 597 
o 

watt - min/in , and b^ = .50 inch. The extent of the 
fusion zone may depend somewhat on the parameter El/Vt 
but it is governed primarily by the thickness t of the 
plate and the geometry of the edge preparation. 

Figures 4 through 10 provide mechanical properties 
at both room temperature and elevated temperatures. Figure 
4 provides typical tension and compression stress strain 
curves for the experimental plate. Figure 5 provides 
the average stress strain curves for varying temperatures 
from room temperature to 700°F based on a half hour 
exposure to the temperature in question. This may or may 
not be a valid approximation for the welding problem, but 
it is the only information currently available. Figure 6 
shows the effect of this temperature on the 0.2% offset 
yield stress. This may be approximated as the cubic 
polynomial given in Table 5. Figure 7 demonstrates the 
effect of the half hour soak temperature on Young's 
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modulus. This curve was generated from Figure 5. Its 
cubic approximation is listed in Table 5. 

Figure 8 shows the room temperature tangent modulus 
at various stress levels. Its straight line approximation 
is 

H' x 10 -3 (ksi) = 0.917a (ksi) + 38.064 

(33) 

The highest stress level shown on this figure is 40.7 ksi. 

+3 

This equates to H' = 0.746 x 10 ksi and may be utilized 

effectively for constant strain hardening. Figure 9 
represents an estimate of the temperature effects on the 
tangent modulus. If one assumes that 



H ' (T°F) « E (T°F) 



(34) 



then 



H » (t°F) 



H' (R.T. ) x E (T°F) 
E (R. T. ) 



(35) 



Masubuchi and Andrews (11) define the constant strain 
hardening parameter for their one-dimensional computer 
program as 

H' = mE (36) 



From the above equations, 



m = H'/E = H' (T°F)/E(T°F) = 0.0742 (37) 



This is the first such value ever determined for their 



Strtu (K*t ) 
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rigor* 7: E*tla*t*d Effect *f T*ap*r*tor* on Yotmg 9 * 

Modulus for 6061|n T6 Condition (31,32). 




<r*i> 
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E««f4 an « Hdlf Hour Ifc p oou r # *t Working 
Tonyrrattiro (32)*. 
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Table 5: Summary of Mechanical and Phyaicai Propertiea for Computer Simulation of 6061 T6 Aluminum Alloy 

2 3 

In the Form P(x) - a ♦ bx ♦ cx ♦ dx * Uae Caution In Reducing the Number of Significant Figures 
Given for the Coefficients at High Temperatures* 



Stress (Ksl) 
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Figure 10 1 Simulation ef Tensile Properties st Various Tssq>eraturea 
after 0.3 Hour st Temp e rature. In this example T°F 
Is 300°P. 



! 
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program. 

Figure 10 also demonstrates a computer algorithm 
for the tensile properties at various temperatures based 
on the half hour soak. To determine this approximate 
value, proceed as follows: 

a. Determine ay(T°F) from Table 5 or Figure 6. 

b. Determine E(T°F) from Table 5 or Figure 7. 

c. Determine H' (T°F) from equation (35) and Table 5 
or Figure 9. 

d. Plot o ( T°F) at 0.2% offset. 

e. Draw line through point d above with slope 
H' (T°F) . 

f. Connect line in e above with line through origin 
with slope E (T°F) . 

5 . Physical Properties 

Figures 11 through 14 show the effect of temperature 
on the coefficient and mean coefficient of linear expansion 
a, density p, thermal conductivity k, and specific heat 
C. The coefficient of linear expansion shown in Figure 11 
is based on empirical equations for high purity aluminum 
and alloy constants provided by reference (35). The alloy 
constant for 6061 TO is 0.990, i.e., 0.990 a (pure 
aluminum) = a (alloy) . The value of this constant for 
heat treated tempers is approximately 0.015 greater. 

With these tempers, application is restricted to temper- 
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atures which do not appreciably exceed those used in the 
final aging treatments (320-350°F). The basic alloy 
constant is limited to temperatures below 600°F. Figure 11 
is extrapolated past 600°F to the melting point. Figure 12 
provides an estimate of density provided by reference 
(36) extrapolated to the melting point. Figure 13 provides 
an estimate of thermal conductivity. Reference (37, 

Volume II, Part 2, pp. 769) provides thermal conductivity 
information of 6000 series aluminum alloys from 300 to 
650°K. A linear relationship was plotted parallel to 
this data passing through the one data point provided by 
reference (38). Reference (37) (Volume I, pp. 11) listed 
the specific heat-temperature relationship for pure aluminum. 
A parallel relationship (39, pp. 180) passing through the 
data point provided by reference (31, Code 3206, pp. 1) 
resulted in the estimated data shown in Figure 14. 

These physical values may be approximated by the 
cubic polynomials given in Table 5 for computer simulation. 
Appendix A lists the computer program for determining 
the coefficients given in Table 5 and for testing these 
coefficients in 10°F increments. 

Table 6 summarizes the physical and mechanical pro- 
perties of 6061 T6 in the form utilized by Masubuchi and 
Andrews (11) in their one-dimensional program. 



Coefficient of Thoraal EKpemeloa (Mlcroecreln/°f) 
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ThtrMl Conductivity (Btu/Hr Ft °T x 10 
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Plgur* 13: (tfiHiM EJ fact of T — p orituro tm Thoraal 
eontocctvtty for 6M1 In T* Cm* It ton. 



Specific |t|t (Mu/Lb^ °F) 
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III: PROCEDURES 



A. Scope of the Research : 

A series of six experiments measuring temperature 
and strain changes during welding was performed. System 
models representing constrained bead-on-plate and butt 
joints were constructed from 6061 aluminum alloy in the 
T6 condition. Welding procedures were utilized which 
offered full penetration with minimum heat input. This 
allowed strain gage locations as close to the weld 
line as possible without exceeding the 400°F maximum 
allowable gage temperature limitation. 

B. Selection o f Parameters : 

The aluminum selected for this investigation was 
picked because of its easy availability and wide use in 
the marine and aeronautical industries. It develops its 
strength from its heat treatment. The 0.250 inch plate 
thickness used in this study was determined by availability, 
past experimental work, ease of handling, and an interest 
in approaching a two-dimensional problem, yet minimizing 
the extent of bending and buckling. The joint design 
for the butt welds used in this experiment (straight) 
was selected for simplicity and ease of preparation. 

Steel backing plates were utilized below the weld and at 
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the clamped ends of the plates to provide support and 
distributed clamping and yet minimize heat conduction 
effects. Each specimen was clamped at its edges to a 
0.500 inch mild steel bed plate which resisted deformation. 

The size of the test plates was set at a nominal 18 
by 30 inches to provide essentially steady state conditions 
at strain gage locations, and to approach the temperature 
distribution of an infinite plate. 

The weld process used was semi-automatic Gas Metal-Arc 
(GMA or MIG). It allows control of weld variables, reduces 
operator error, and fosters repeatability. No preheating 
was necessary. 

C. Strain Measurement by Electric Resistance Strain Gages: 

The fundamental concept of strain gage operation 
is that certain conductors exhibit a change in electrical 
resistance with a change in strain. Gages designed 
according to this principle are attached to test materials 
whose strains are then monitored by measuring resistance 
variations across the gage. In the case of welding 
thermal strains, the observed resistance change, AR» is 
made up of : 



AR = AR 1 Ce e ) + AR 2 ( £p) + AR 3 (aT) + AR 4 (T) 
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where 






Ar 3 (aT) 

ar 4 (T) 



the resistance change corresponding to 

elastic mechanical strain, e , from which 

e 

stresses can be computed. 

the resistance change corresponding to 

plastic mechanical strain, e , if it exists. 

P 

the resistance change corresponding to 
temperature induced thermal strain, aT. 
the resistance change caused by thermo- 
electric effects in the gage itself. 



While it is not presently possible to discriminate 
between the two mechanical strains, e and e , Ar (aT) 

© P J 

and AR^ (T) can be separated out by physical temperature 
compensation of the strain gage material coupled with 
empirical calculation. For this investigation, a T for the 
gage was approximately 13.0 microstrain/°F . The AR 4 (T) 
and the difference between gage a T and base plate a^ 
will then generate an error. To compute this error, 
a test gage of the type and lot used in the experiment 
is mounted on a small sample of 2024 aluminum alloy by 
the gage manufacturer. The sample is then heated at 
equilibrium until a curve of "apparent strain" vs. temper- 
ature is obtained to the operating temperature range. 

The gage readings recorded in the weld experiments can 
then be corrected by subtracting out the apparent strain 
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value corresponding to temperatures observed at the gage 
location. Figures 16 and 17 show this factory provided 
information for the individual gages and rosettes 
respectively. In each case, the curves have been 
reduced by the factory to fourth order polynomials yield- 
ing microstrain for ease in computer calculation. 

Gage: 

EAP = -69.36 + 2.07T - 2.02 x 10~ 2 T 2 + 5.39 x 10 _5 T 3 

- 3.27 x 10 _8 T 4 (38) 



Rosettes : 

EAP = -60.68 + 2.22T - 2.40 x 10~ 2 T 2 + 7.05 x 10 _5 T 3 



-8 4 

- 5.59 x 10 T 



(39) 

Since 



6061 T6 has a greater coefficient of thermal expansion 
than the 2024 T4 test plate, a correction to these poly- 
nomials has been made. This correction was calculated 
as follows: 

In general, 



Thermal Strain = a(T) dT 

Jt ! 

From reference (35) , 

“ 6061 = c 6061 “ pure A1 - = .99 a pure Al. 



“ 2024 “ C 2024 “ P ure A1 



= .97 a pure Al. 



a 



a pure Al 

T„ 



6061 

.99 



AEAP = 



(c 6061 ~ c 2024> a pure Al. dT 



(40) 

(41) 

(42) 

(43) 

(44) 
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AEAP = 




(45) 



Integrating the polynomial given in Table 5 for a 



6061 <T) 



0 <_ T _< 400°F, between 77°F and T°F yields the following 
apparent strain correction for both gages and rosettes: 



Equation (46) must be added to equations (38) and (39) 
to obtain true apparent strain corrections. The true 
apparent strain is then subtracted from the actual strain 
to provide true mechanical strains. 

Strain gages manufactured in the form of rosettes 
provide as many as three independent readings in three 
directions at a single location. This is enough to 
completely describe the two-dimensional strain state 
at that location (41, 42, 43) and is discussed in the 
next chapter. 

D. Description of Apparatus: 

1. Test Plates: A description of the test 

specimens is summarized in Figures 18 through 22 which 
provide dimensions and gage locations. Devices Q) through 
(?2)were dynamically measured during welding. Gages 

through were statically measured before and after 



AEAP = -19.4979 + .2463T + .9561 x 10 4 T 2 
+ .5303 x 10~ 10 T 4 



.7879 x 10 -7 T 3 



(46) 




O Thermocouple 



SCraln Gage 



Strain Gage Rosette 
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Figure 18: Test No. 1 and Test No. 2 (Butt and Bead-on-Plate 

Welds Respectively) Measuring Temperature Distribution. 



7T 



9 * 






Figure 19: Test No. 3 (Bead-on-Plate Weld) Measuring Temperature 

Distribution and Strains. 



© Thermocouple ea. Strain Gage <^Q<7 strain Gage Rosette 

f~i Extensiometer Tab — Tack Weld f 





Figure 21: Test No. 5 (Butt Weld) Measuring Temperature 

Distribution. Strains, and Extension of Root Gap. 
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o Thermocouple «=» Strain Gage ^0^ Strain Gage Rosette 



□ Extensiometer Tab Tack Weld 




Figure 22: Test No. 6 (Butt Weld) Measuring Temperature 

Distribution. Strains, and Extension of Root Gap. 
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welding. The arrangement of constraining clamps is 
shown in Figure 23. 



2. Sensors and Instrumentation ; 
a. Strain Gages. Two types of gages were used 
in this investigation: SR-4 foil 45° rosettes and 

SR-4 foil single element gages. Gage properties 



were as follows: 
Gage 

Designation 
Manufacturer 
Grid Length (in.) 

Grid Width (in.) 
Overall Length (in.) 
Overall Width (in. ) 
Temperature Range 
Resistance 
Gage Factor 
Cement 

Protective Covering 



SR-4 Rosette 
FAER-18RB-12S13-ET 
BLH Electronics 
3/16 
0.90 
0.280 
0.540 

-50 to +400°F 

120 . 

2.03 
EPY-600 
BLH Barrier-C 



SR-4 Gage 
FAE-25-12S13L 
BLH Electronics 
1/4 
0.13 
0.35 
0.13 

-50 to +400°F 
120 . 

2.07 
EPY-600 
BLH Barrier-C 



b. Temperature Sensors: All temperature sensors 

used were Chrome 1/Alumel thermocouples made from 
Leads and Northrup No. 28 wires. Each thermocouple 
was spot-welded onto the test specimen and protected 
by No. 33 Sauereisen Sealing Cement. Each was 
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Figure 23: Constraining Equipment 
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located 0-536 inch (1 second of arc travel) ahead 
of the strain gage mid -axis. 

c. Instrumentation. Strain gages were connected 
into a Potentiometric Circuit (Half-Wheatstone 
Bridge) , balanced and calibrated as indicated 
schematically in Figure 24. Thermocouples were 
referenced to a 32°F ice-bath and calibrated as 
indicated in Figure 25. Both circuits were fed into 
a Honeywell continuous-recording, 12-channel 
Visicorder. When the raw data was actually read 
off the recorder tape, some traces were delayed or 
advanced with respect to others. This was done 

to correct for the finite difference in position 
along the weld line of the thermocouples (1 second) 
and strain rosette elements (+ 1/3 second) . 

This produces the effect of a simultaneous 
strain and temperature reading at the strain gage 
location. Timing was accomplished by means of an 
electric stop watch as well as a timer integral 
to the Visicorder. Static strain measurements 
were taken with a Budd Instrument Division Model 
P-350 strain indicator. 

d. Extensiometer. Figure 26 shows a sketch of 

the extensiometer utilized in Tests 5 and 6. Strain 
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Battery 



Strain Gage Instrumentation Circuit (7.) 



Wheelco Potentiometer 
Model No. 320P 




Junction 



Figure 25: Thermocouple Instrumentation Circuit (7) 
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Figure 26: Ex ttmloMttr* The Adjusting 

Screws art Set So That the Spring Holds 
Extensiometer Firmly In Place. 




Figure 27: 



CxttnsloMCtr Calibration Curve. 0.000 Inches Represents 
an Extensiometer Gap of 3.00 Inches. 
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gages were located on both sides of a spring steel 
ribbon and connected into adjacent legs of the 
resistance network circuit described in paragraph 
c. above. The net effect of such a set-up is to 
double the bending sensitivity and null the stretching 
sensitivity. Figure 27 provides the calibration 
curve for the Visacorder referenced to actual 
dimensional changes in inches. The extensiometer and 
its mounting tabs may be seen in the photographs 
in Figures 28 through 31. It was pulled back away 
from the torch to prevent damage when the arc past 
the mounting tabs. Once the torch was safely clear, 
it was set back in place to observe weld shrinkage. 

It is noted that the mounting tabs were cemented in 
place with EPY 600 cement vice bolting or welding 
to minimize local straining of the base plate. 

3. Welding Equipment and Conditions : Welding conditions 

are summarized in Tables 7 and 8. The welding machine 
utilized was manufactured by the Linde Division of Union 
Carbide Corporation and consisted of a type HW-16 GMA 
torch, a SVI-300 power supply, and associated governor, 
carriage and wire feed mechanisms. Travel speed, arc 
voltage, amperage and arc length were present to the same 
value for all experiments. Some fluctuations of amperage 
did exist when the geometry of the weld changed from 




Figure 28 2 Overview of Exper imental Equipment Showing Instrumentation and Recorder 
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Figure 30! Plate Instrumented and Clamped Under Welding Torch 
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Figure 31; Sample Plated Instrumented and Welded 
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bead-on-plate to butt. Wire feed is a functin of the 
other variables and was maintained by the machine auto- 
matically* Figures 29 through 31 are photographs of 
the actual equipment set-up. 

E. Experimental Procedures: 

The experimental operation is shown schematically in 
Figure 32. The test plate was instrumented and clamped 
into place. The welding machine was lined up with the 
joint and positioned over the run-off tab at the left 
end of the backing plate. Welding speed, arc voltage and 
amperage were pre-set. The visicorder was actuated and 
an arc was struck on the backing plate. As the welding 
torch crossed from starting tab to plate, the timer was 
started. The Visicorder output was marked when the arc 
passed the center strain gage location. When the welding 
bead reached the run-off tab at the right end of the 
plate, the arc was extinguished and the plate allowed to 
cool. Thr recorder continued to monitor the gages until 
the plate cooled to ambient temperature. After cooling, 
the plate was released from its clamp, and the new strains 
recorded either with the recorder or the Budd instrument ^ 




Figure 32: Schematic of Apparatus and Procedure. (7) 
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IV: DATA DEDUCTION 



A. General: 

At any point on a free (unloaded) surface of a 
solid, it is necessary to know three independent quantities 
in order to specify the state of stress completely. These 
quantities are the magnitudes of two principal stresses, 
o^, and and their directions, (J) or (<|> + 90°), with 

respect to some reference. For isotropic elastic materials, 
these values can be calculated from strains measured on 
the surface at the point in question, and since three 
independent quantities are to be determined, it will be 
necessary to make three independent measurements of strain 
(42, pp. G2-03) . 

B. Strain Calculations: 

Equation 47 describes the relationship for general 
strains for a three-element rosette shown in Figure 33. 




Figure 33 
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In the special case of a three-element 45° rectangular 



rosette shown in Figure 34, 
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Equation (47) reduces to 
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Solving Equation (48) for e , e , and y , 

x y xy 
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(49) 



The Mohr's Circle for this special case is shown in Figure 35, 
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Figure 35: Mohr's Circle for the Rectangular Rosette With 

Three Observations of Strain (42 t pp» G2-12). 



Principal strains are 
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where 4> is defined as the angle of one of the principal 
axes with respect to the axis of reference. Mathematically, 



Tan 2<f> 




(52) 



In the literature, there is considerable confusion 
concerning 4> and which principal axis it defines. 

Murray (42, pp. G2-11, 15) offers a consistent method 
for dealing with this problem. If 4>^is defined as the 
angle measured (positive in the anti-clockwise direction) 
from the positive OA axis of the strain rosette to the 
positive 01 axis which corresponds to the direction of 
e .. , the algebraically largest principal strain, then: 
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Rule 3- When 2 — - (55) 

and (a) e > e then e = £ and 4> = 0° (55a) 

a. C a. -L -L 

or (b) e < e c , then ^ = e 2 and 4> 1 = + 90° (55b) 
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Rules 1 through 3 are defined with axis OA as the 
reference. To convert to the situation consistent with 
this investigation, axis OA is +45° (anti-clockwise) from 
the x-axis. 

B 




4> x = $1 + 45 ° (56) 

The geometry utilized in this investigation reduces 
to that of the mathematical model if the axes are kept 
consistent with Figure 36. 

C. Stress Calculations: 

This investigation yields mechanical strains. The 
automatic temperature compensating feature of the strain 
gage coupled with the apparent strain correction removes 
the thermal strain effect. Stresses must be determined 
by total strain. 

e mech + e thermal ~ e total (57) 



where 
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For general two-dimensional stress calculations. 
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For this particular configuration, 
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v is Poisson's ratio. Nominally, v = .333 for aluminum. 
v will probably change with temperature, but should not 
exceed a maximum value of .500 in the weld puddle because 
of volume consideration when melted (at .500, there will 
be no volume change) . 

Stress calculations are dependent upon E(T) and a (T) 
assuming v is relatively independent of temperature. These 



